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ABSTRACT The study of passive scalar transport in a turbulent 
velocity field leads naturally to the notion of generalized flows which 
are families of probability distributions on the space of solutions to 
the associated ODEs which no longer satisfy the uniqueness theo- 
rem for ODEs. Two most natural regularizations of this problem, 
namely the regularization via adding small molecular diffusion and 
the regularization via smoothing out the velocity field are consid- 
ered. White-in-time random velocity fields are used as an example to 
examine the variety of phenomena that take place when the velocity 
field is not spatially regular. Three different regimes characterized by 
their degrees of compressibility are isolated in the parameter space. 
In the regime of intermediate compressibility, the two different regu- 
larizations give rise to two different scaling behavior for the structure 
functions of the passive scalar. Physically this means that the scal- 
ing depends on Prandtl number. In the other two regimes the two 
different regularizations give rise to the same generalized flows even 
though the sense of convergence can be very different. The "one 
force, one solution" principle is established for the scalar fleld in the 
weakly compressible regime, and for the difference of the scalar in the 
strongly compressible regime which is the regime of inverse cascade. 
Existence and uniqueness of an invariant measure is also proved in 
these regimes when the transport equation is suitably forced. Fi- 
nally incomplete self-similarity in the sense of Barenblatt-Chorin is 
established. 



Introduction 

Recent efforts to understand the fundamental physics of hydrodynamic turbu- 
lence have concentrated on the explanation of the observed violations of Kol- 
mogorov's scaling. These violations reflect the occurrence of large fluctuations in 
the velocity field on the small scales, a phenomenon referred to as intermittency. 
Some progress in the understanding of intermittency has been achieved recently 
through the study of simple model problems that include Burgers equation 
[1,0 and the passive advcction of a scalar by a velocity field of known statistics 
[3,^ ^, ^. This paper is a summary of the many interesting mathematical 
issues that arise in the problem of passive scalar advection together with our 
understanding of these issues. We put some of our results in the perspective 
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of a new phenomenological model proposed recently by Barenblatt and Chorin 
H using the formalism of incomplete self-similarity. 

Generalized Flows 

Consider the transport equation for the scalar field 6'^{x, t) in M.'^: 

— + {u{x,t)-V)d'' = KM\ (1) 

We will be interested in 0" in the limit as k — > 0. It is known from classical 
results that if u is Lipschitz continuous in x, then as k — > 0, 6'^ converges to 6, 
the solution of 

— + {u{x,t)-V)e^O. (2) 

Furthermore, if we define {ipg f {x)} as the flow generated by the velocity field 
u, satisfying the ordinary differential equations (ODEs) 

— — =u{ips,t{x),t), ipsA^)=x, (3) 

for s < t, then the solution of the transport equation in ^ for the initial condition 
6'^{x, 0) = 6*0(2:) is given by 

0{x,t) = eo{^o,li^))^MvtA^))- (4) 

This classical scenario breaks down when u fails to be Lipschitz continuous in 
x, which is precisely the case for fully developed turbulent velocity fields. In 
this case Kolmogorov's theory of turbulent flows suggests that u is only Holder 
continuous with an exponent roughly equal to ^ for c? = 3. In such situations 
the solution of the ODEs in || may fail to be unique |^ , and we then have to 
consider probability distributions on the set of solutions in order to solve the 
transport equation in bl This is the essence of the notion of generalized flows 
proposed by Brenier jl^, |l^ (see also [l3|]). 

There are two ways to think about probability distributions on the solu- 
tions of the ODEs in |. We can either think of it as probability measures on 
the path-space (functions of t) supported by paths which are solutions of or 
we can think of it as transition probability at time t if the starting position 
at time s is x. In the classical situation when u is Lipschitz continuous, this 
transition probability degenerates to a point mass centered at the unique so- 
lution of ^ When Lipschitz condition fails, this transition probability may be 
non-degenerate and the system in |^ is intrinsically stochastic. 

There is a parallel story for the case when w is a white-in-time random 
process defined on a probability space {^1,^,1^). We will denote the elements 
in 17 by w and indicate the dependence on realization of the random velocity 
field by a super- or a subscript lo. In connection with the transport equation in 
H, it is most natural to consider the stochastic ODEs 

dftti^) = '"(¥'^(2;),*)^^, flsix) = X, (5) 
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in Stratonovich sense. In this case, it is shown that if the local characteristic 
of u is spatially twice continuously differentiable, then the system in || has a 
unique solution. Such conditions are not satisfied by typical turbulent velocity 
fields on the scale of interest. When the regularity condition on u fails, there 
are at least two natural ways to regularize || or ||. The first is to add diffusion: 

d<f (a:) = (x), t)dt + ^/2^dP{t), (6) 

and consider the limit as k ^ 0. We will call this the K-limit. The second is to 
smooth out the velocity field. Let tpi, be defined as i^eix) = e~'''ip{x/e), where 
ip is a standard moUifier: > 0, J^j^ipdx — 1, decays fast at infinity. Let 
■u^ — ui^ -tp^ and consider 

(7) 

in the limit as £ ^ 0. We will call this the e-limit. Physically k plays the role of 
molecular diffusivity, e can be thought of as a crude model of the viscous cut- 
off scale. The K-limit corresponds to the situation when the Prandtl number, 
defined here as the ratio of e and k, tends to zero, Pr — + 0, whereas the e-limit 
corresponds to the situation when the Prandtl number diverges, Pr —^ oo. The 
following questions naturally arise: 

(Ql) How do the flows and the passive scalar behave statistically in the k- and 
£-limits? 

(Q2) Does there exist a unique statistical steady state when the transport equa- 
tion in |l] is suitably forced? 

(Q3) What are the statistical and geometrical properties of solutions in the 
statistical steady state? 

Below we address these questions on a specific model introduced by Kraichnan 
[I- 

Before proceeding further, we relate the regularized flows in to the 
solutions of the transport equations. Consider the K-regularization first. It is 
convenient to introduce the backward transition probability 

g^{x,t\dy,s)^Ep6{y-^ts"{x))dy, s < t, (8) 

where the expectation is taken with respect to P{t), and ip'^]^{x) is the fiow 
inverse to ip^'^{x) defined in |^ (i.e. ip^'f{x) is the forward flow and (p^'^{x) is 
the backward flow). The action of g'^ generates a semi-group of transformation 

/ ^{y)gZ{x,t\dy,s), (9) 

JR<^ 

for all test functions 0^{x,t) = S'^'^^{x) solve the transport equation in |^ 
for the initial condition 0^(x, s) = il>{x). Similarly, for the flow in ^ define 

= s<t. (10) 
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6^{x,t) = S^fxpix) solves the transport equation 

^ + K(x,t).V)r = 0, (11) 

with initial condition s) — ip{x). Similar definitions can be given for forward 
fiows but we will restrict attention to the backward ones since we are primarily 
interested in scalar transport. The results given below generalize trivially to 
forward flows. 



Kraichnan Model 

In |l5) Kraichnan introduced one of the simplest model of passive scalar by con- 
sidering the advection by a Gaussian, spatially non-smooth and white-in-time 
velocity field. The fact that white-in-time velocity fields may exhibit intermit- 
tency was first recognized by Majda [l6| pTf . Definitive work on Kraichnan 
model has been done afterwards in [|[ 

We will consider a generalization of Kraichnan model introduced in | [T^ (see 
also ||l^). The velocity u is assumed to be a statistically homogeneous, isotropic 
and stationary Gaussian field with mean zero and covariance 

'E,Ua{x,t)uf3{y, s) = {CaSap - Caf3{x - y))S{t - s). (12) 



We assume that u has a correlation length £q, i.e. the covariance in ^ decays 
fast for |x — ?/| > £o. Consequently Ca/3{x) — > CqiSq/j as |a;|/£o oo. On the 
other hand, we will be mainly interested in small scale phenomena for which 
\x\ ^ £o. In this range, we take Cai3{x) — dap{x) + 0{\x\^ /£q) with 

d^p{x) = Adlp{x) + Bdip{x), (13) 

and 



dL{x) = I? (d + e - i)5o.p - 1 bl« 



(14) 



D is a parameter with dimension [Icngth]^^^ [timc]^^. The dimensionless pa- 
rameters A and B measure the divergence and rotation of the field u. ^ = 
corresponds to incompressible fields with V • m = 0. B — corresponds to irro- 
tational fields with V x u = 0. The parameter f measures the spatial regularity 
of u. For ^ S (0,2), the local characteristic of u fails to be twice differentiable 
and this fact has important consequences on both the transport equation in |^ 
and the systems of ODEs in || or 

Existing physics literature concentrates on the K-limit for Kraichnan model. 
Let = A + {d- 1)B, = A,V = C^/S^. V € [0, 1] is a measure of the 
degree of compressibility of u. The pioneering work of Gaw§dzki and Vergassola 
[ pis) (see also |l^) identifies two different regimes for the K-limit: 



4 



1. The strongly compressible regime when V > d/^^. In this regime 
converges to a flow of maps, i.e. there exists a two-parameter family of maps 
{<p"j,(x)} such that 

g::ix,t\dy,s)^6{y-^^Jx))dy. (15) 

Moreover particles have finite probability to coalesce under the flow of {Lpfg{x)}. 
In other words the flow is not invertible. 

2. When V < d/^^, converges to a "generalized stochastic flow" 

g^{x,t\dy,s) ^ g^{x,t\dy,s), (16) 

and the limit is a nontrivial probability distribution in y. This means that 
the image of a particle under the flow defined by the velocity field u is non- 
unique and has a non-trivial distribution. In other words, particle trajectories 
branch. 

The same classification of the fiows was obtained by Le Jan and Raimond 
using Wiener chaos expansion without explicit reference to the k limit. In 
contrast, our primary motivation is to study the limit of physical regularizations. 

The following result answers the question Ql and also points out that there 
are three different regimes if both the k- and the e-limits are considered. 

Theorem 1 In the strongly compressible regime when V > d/^^, there exists a 
two-parameter family of random maps {ip'^g{x)}, such that for all smooth test 
functions tp and for all (s, i, x), s < t, 

E(5-;Xx)-^«,(x)))'^0, (17) 

as K — > 0, and 

E{i.{^l':{x))-i;{^-Jx))y^O, (18) 

as e — )■ 0. Moreover, the limiting flow {^Pt si^)} coalesces in the sense that for 
almost all (t,x,y), x ^ y, we can define a time t such that — oo < t <t a.s. 
and 

Vlsi^) ^ VUV) for s<T. (19) 

In the weakly compressible regime when V < (d + ^ — 2)/2^, there exists a 
random family of generalized flows guj{x, t\dy, s), such that for all test function 

S^,Mx)= f ij{y)gUx,t\dy,s), (20) 

satisfies 

E{Stf^P{x)-S-M^)y^O, (21) 
as K — !■ for all (s, t, x), s < t, and 

e([ rjix){^l,ipl':{x))-Sr.M^))dxY ^0, (22) 
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as e — *■ for all {s , t) , s <t, and for all test functions r] . Moreover, g^(x,t\dy, s) 
is non- degenerate in the sense that 

Sr^,,p^{x)-{S'^,M^)y>0 a.s. (23) 

In the intermediate regime when (rf + ^ — 2)/2^ < V < d/C^, there exists a 
random family of generalized flows gi^{x, t\dy^ s), such that for all test function 
V' and for all {s,t,x), s <t, 

E{srf^P{x)-SZ,^P{x)f ^0 (24) 

as K — !■ 0. In the e -limit, the flows ip^f{x) converges in the sense of distribu- 
tions, i.e. there exists a family of probability densities 

{Gn{xi, . . .,Xn,t\yi, . . .,yn,s)dyi ■ --dyn}, (25) 

n ~ 1,2, . . ., such that 



(26) 



'E''iPi'Pt,fixi),- ■ ■ ,ipt '^%Xn)) ^ J ipiyi,- ■ ■ ,yn) 
xGnixi, ■ ■ ■ ,Xn,t\yi, . . . ,yn, s)dyi ■■■dyn, 



as e for any continuous function ip with compact support. Furthermore, 
the limit coalesces in the sense that 

G2{xi,X2,t\yi,y2,s) ^ G2{xi,X2,t\yi,y2,s) 

(27) 

+A{yi,xi,X2,t, s)S{yi - 2/2), 

with yl > when t > s. Here G2 is the absolutely continuous part of G2 with 
respect to the Lebesgue measure. Similar statements hold for the other Gn 's. In 
particular, the {Gn} 's differ from the moments of the n-limit g^^ defined in |^. 

Rephrasing the content of this resuh, we have strong convergence to a family 
of flow maps in the strongly compressible regime for both the K-limit and the 
e-limit. In the weakly compressible regime, we have strong convergence to a 
family of generalized flows for the K-limit, but weak convergence to the same 
limit for the e-regularization. In fact, using the terminology of Young measures 
[ pO| , the limiting generalized flow {g^ix, t\dy, s)} is nothing but the Young mea- 
sure for the sequence of flow maps {i^a^f (x)}. Finally, in contrast to what is 
observed in the other two regimes, the e-limit and K-limit are not the same in 
the intermediate regime. As we will see below, the structure functions of the 
passive scalar field scale differently in the two limits. 

From Theorem ^, it is natural to define the solution of the transport equation 
in 1^ for the initial condition 0^{x, s) = 0o{x) as 

e^{x,t) ^ S'^Mx) = I eo{y)g^{x,t\dy,s), (28) 
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for the weakly compressible and the intermediate regimes in the K-limit (non- 
degenerate cases), and as 



eUx,t) = OoiifiLix)), (29) 

for the strongly compressible regime. In the intermediate regime in the e-limit, 
it makes sense to look at the limiting moments of 6^ {x, t) since we have as e ^ 



E(0^(xi,i)---C(^n,i))^ / Ooiyi)---0oiyn) 

xGn{xi, ■ ■ ■,Xn,t\yi, . . .,yn,s)dyi ■ ■ ■ dy^i. 

It should be noted that when is non-degenerate, there exists an anomalous 
dissipation mechanism for the scalar, whereas no such anomalous dissipation is 
present in the coalescence cases The presence of anomalous dissipation 

is the primary reason why the transport equation in |^ has a statistical steady 
state (invariant measure) if it is appropriately forced, as we will show later. 

Details of the proof of Theorem |l| are given in [|l] . Crucial to the proof is 
the study of P{p\r, s) defined through e-regularization as 



T^{r)P{p\r,s~t)dr^\iin^rj{\^'t::{y)~vts{^^)\), (31) 

£—♦0 ' 



where ry is a test function, and similarly through K-regularization. Here p = 
\y — z\ and s < t. P{p\r,s) can be thought of as the probability density that 
two particles have distance r at time s < t if their final distance is p at time t. 
For Kraichnan model, P satisfies the backward equation 
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for the final condition Yivcis^o^ P[p\r^ s) = 6{r — p), and with a(r), b{r) such 
that 

air)=D{S' +^C'y + Oir'/il), 
h{r) = D{{d - 1 + - eC2)r«-i + 0{r/tl). 
For r ^ (.Q, a(r) tends to Co, b{r) to Co((i — l)/?^, and the equation in 32| 
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reduces to a diffusion equation with constant coefficient. The equation in 32 
singular at r = 0. The proof of Theorem ^ is essentially reduced to the study 
of this singular diffusion equation. This is also the main step for which the 
white-in-time nature of the velocity field is crucial. 



Structure Functions 



We now study some consequences of Theorem |l] for the passive scalar 9ui defined 
in E8^ or 29 . We note that the scaling of the second-order structure function is 



the same for the k- and the e-limits in the strongly and the weakly compressible 
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cases [£2j, but it differs in tfie intermediate regime as a result of tfie difference 
between tlie limits in ^ and |2^. For simplicity of presentation, we assume that 
^0 is isotropic and Gaussian. Denote (n € N) 

S2n{\x - y\,t) = E{0Ux, t) - e^{y,t)f^, (34) 

or S2n{\x~ylt)^Yimn6l{x,t)-ei{y,t)f^, (34') 

£—►0 

in the intermediate regime in the e-limit. In the strongly compressible case, we 
have for both the n- and the e-limits 

52(r,i) = 0(/), (35) 

with 

In the weakly compressible case, we have for both the k- and the e-limits 

^2(r,i)=0(r2-«). (37) 

In the intermediate regime, the limits differ, and the K-limit scales as in 

can be derived 



whereas the e-limit scales as in The equations in |35| and |37 
upon expressing S2 in terms of P; the details are given in Ref . 1 2l| . 

It is interesting to discuss the higher order structure functions both in the 
non-degenerate and in the coalescence cases in ^ and ^ since their scalings 
highlight very different behavior of the scalar. We consider first the coalescence 
cases which are simpler. In these cases, because of the absence of dissipative 
anomaly, all higher order structure functions can again be expressed in terms 
of P, and it can be shown that 

52„(r,i) =0(/), (38) 

with C given byjs^ for all n > 2. In fact, coalescence implies that the temperature 
field Ouj tends to become flat except possibly on a zero-measure set where it 
presents shock-like discontinuities. Such a situation with two kinds of spatial 
structures for 9^ is usually refered to as bi-fractal, and, in simple cases, one may 
identify ^ with the codimension of the set supporting the discontinuities of O^j 

[ilil- 

The non-degenerate cases are more complicated. In these cases, one expects 
that 9^^ presents a spatial behavior much richer than in the coalescence cases, 
with all kinds of scalings present. This is the multi- fractal situation for which 
the higher order structure functions behave as 

52„(r,t) = 0(r«-), (39) 

with C,2n < "-(2 — £,) for 2n > 2. The actual value of the Cn's cannot be ob- 
tained by dimensional analysis, and one has to resort to various sophisticated 
perturbation techniques (see |||, |^, |[ ^). We will consider again the scaling of 
the structure functions at statistical steady state in the section on incomplete 
self-similarity. 
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One Force, One Solution Principle for Tempera- 
ture 



We now turn to question Q3 and consider the existence of a statistical steady 
state for the transport equation with appropriate forcing. We restrict atten- 
tion to the non-degenerate cases which include the weakly compressible regime 
and the intermediate regime in the K-limit. Indeed, in these regimes the non- 
degeneracy of gui{x, t\dy, s) as a probability distribution in y implies dissipation 
of energy or, phrased differently, decay in memory in the semi-group St^s gener- 
ated by {Qlj}. We show that the anomalous dissipation is strong enough in order 
that the forced transport equation has a unique invariant measure for both the 
weakly compressible regime and the intermediate regime in the K-limit. This 
result, however, depends on the finiteness of Ho. In limit as ~* oo an invariant 
measure exists only for the weakly compressible regime. 
We will consider (compare with ^) 

pin 

— + (u{x,t)-\/)0^b{x,t). (40) 

where 6 is a white-noise forcing such that 

i:b{x,t)b{y,s)^B{\x-y\)S{t-s). (41) 

B{r) is assumed to be smooth and rapidly decaying to zero for r ^ L; L will 
be referred to as the forcing scale. The solution of ^ for the initial condition 
0^{x,s) — Oo{x) is understood as 



9Ux,t) = Srjoix) + Sl,b{x,T)dT. (42) 

J s 

Define the product probability space (ri„ x fib, x J-h,Vu x Vb), and the shift 
operator TT-Lu{t) = Lo{t + r), with cu = {lOu, oJb)- We have 

Theorem 2 (One force— one solution I) For d > 2, in the weakly compress- 
ible regime and in the intermediate regime in the n-limit, for almost alluj, there 
exists a unique solution of^^ defined on M'' x (—00,00). This solution can be 
expressed as 

0*Jx,t) = / SrMx,s)ds. (43) 



Furthermore the map ui ^ O'^ satisfies the invariance property 

0*TU^,t)^eUx,t + r). (44) 

Theorem 1^ is the "one force, one solution" principle articulated in Because 
of the invariance property^, the map in pleads to a natural invariant measure. 
As a consequence we have 
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Corollary 3 For d > 2, in the weakly compressible regime and in the intermedi- 
ate regime in the n-limit, there exists a unique invariant measure on Lf^^(M.'^ x 17) 
for the dynamics defined by 



The connection between the map 43 and the invariant measure, together with 
uniqueness, is explained in |2^. The restriction on the dimensionality in Theo- 
rem || arises because the velocity field has finite correlation length Iq: Theorem]^ 
is changed into Theorem ^ below in the limit as -^o — oo which can be considered 
after appropriate redefinition of the velocity field. 

We sketch the proof of Theorem |[ Basically, it amounts to verifying that 
the dissipation in the system is strong enough in the sense that 



E 



T2 

Ti 



/ 5,%,6(a;,s)ds)'^0, (45) 



as Ti, Ti — s- — CX3 for fixed x and t. The average in ^ is given explicitly by 

/ / B{p)P{^\p,s)drds, (46) 

where P satisfies ^ The convergence of the integral in ^ depends on the rate of 
decay in \s\ of P(0|p, s). The latter can be estimated by studying the equation in 
which yields F(0|p,s) - Cp°'\s\-'^/^ with a = {d^l-^{^+l)V)/ {l+^V) 
for \s\ large and p <^ £q. Hence, the integral in s in ^ tends to zero as Ti, 
T2 — > — c» if d > 2. It follows that the invariant measure in ^ exists provided 
d>2. 

We now ask what happens if we let £0 —^ 00 in order to emphasizes the 
effect of the inertial range of the velocity? This question, however, has to be 
considered carefully because the velocity field with the covariance in ^ diverges 
as £0 ~^ 00. The right way to proceed is to consider an alternative velocity u, 
taken to be Gaussian, white-in-time, but non-homogeneous, with covariance 

Evaix,t)v(}{y,s) 

= {cafsix - a) -\- Ca0{a- y) - Ca/six - y))S{t - s). 

For finite £0, one has v{x,t) = u{x,t) — u{a,t), where a is arbitrary but fixed. 
However, v makes sense in the limit as £q 00. Denote by 'di^{x,t) the tem- 
perature field advected by v, i.e. the solution of the transport equation ^ with 
u replaced by v: 

^-^{v{x,t)-V)i)^b{x,t). (48) 

Restricting to zero initial condition, it follows from the homogeneity of the 
forcing that the single-time moments of 9^ and -duj coincide for finite £0, but in 
contrast to 9i^, d^j makes sense as £q ^ 00. Thus, is a natural process to 
study the limit as ^ 00, and we have 
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Theorem 4 (One force— one solution II) In the limit as £o ^ oo in the 

weakly compressible regime, for almost all uj, there exists a unique solution of 
defined on M."^ X (— oo,c»). This solution can be expressed as 

rjx,t)^ f SfJix,s)ds, (49) 

where S^^ is the semi-group for the generalized flow associated with the velocity 
defined in ^ in the limit as £q ~* oo. Furthermore the map lj -d*^ satisfies 
the invariance property 

rT^Jx,t)=d*Jx,t + T). (50) 
As a direct result we also have 

Corollary 5 In the limit as io ^ oo, in the weakly compressible regime there 
exists a unique invariant measure on Lf (M*^ x il) for the dynamics defined by 

Notice that, as £a oo, the anomalous dissipation is not strong enough in the 
intermediate regime in the K-limit, for which no statistical steady state with 
finite energy exists. 

The proof of Theorem ^ proceeds as the one for Theorem but the estimate 
for P in |§ changes as P(0|p, s) ~ Cp"|s|-("+i)/(2-C) with a = (d - 1 - ^(C + 
1)P)/{1 + ^V) for \s\ large and p <C £q- It follows that the integral in s in ^6] 
converges as Ti, r2 — > — oo in the weakly compressible regime only. 

One Force, One Solution Principle for the Tem- 
perature Difference 

Since no anomalous dissipation is present in the coalescence cases, i.e the stron- 
gly compressible regime and the intermediate regime in the e-limit, no invariant 
measure for the temperature field exists in these regimes. It makes sense, how- 
ever, to ask about the existence of an invariant measure for the temperature 
difference, i.e. to consider 

Se^{x,y,t)^ I StMx,s)^b{y,s))ds, (51) 

JT 

in the limit as T ^ — oo. When 0* exists, one has 

50*(a;,y,t)= hm y, t) = 0:^(x, t) - 0* (y, t), (52) 

1 — > — oo 

but it is conceivable that (56'* exists in the coalescence cases even though 0* 
is not defined. The reason is that coalescence of the generalized flow implies 
that the temperature field flattens with time, which is a dissipation mechanism 
as far as the temperature difference is concerned. Of course, this effect has 
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to overcome the fluctuations produced by the forcing, and the existence of an 
invariant measure such as ^ will depend on how fast particles coalesce under 
the flow, which happens only in the limit as io ^ oo (i.e. for the alternate 



velocity deflned in 47) as we show now. 

For finite Iq, if we were to consider two particles separated by much more 
than the correlation length Iq, the dynamics of their distance under the flow 
is governed by the equation in ^ for r £o, i-e. by a diffusion equation with 
constant diffusion coefficient on the scale of interest. It follows that no tendency 
of coalescence is observed before the distance becomes smaller than io, which, 
as shown below, does not happen fast enough in order to overcome the the 
fluctuations produced by the forcing. In other words. 

Lemma 6 In the coalescence cases, for finite £o, there is no invariant measure 
with finite energy for the temperature difference. 

Consider now the limit as £o oo, and let 5'di^{x,y,i) — z?^(a;,t) — ^^{yjt) 
where solves the equation in The temperature difference d-d^ satisfies 
the transport equation 

dS-d 

— + (vix, t)-V, + v{y, t) ■ Vy)Sd = b{x, t) - b{y, t). (53) 

We have 

Theorem 7 (One force— one solution III) In the limit as Iq oo, for al- 
most alluj, in the strongly and the weakly compressible regimes, as well as in the 
intermediate regime if the flow is non-degenerate, there exists a unigue solution 



of 5S defined on R'^ x (— oo, oo). This solution can be expressed as 

Srjx,y,t)^ f Sr,M^,s)~b{y,s))ds, (54) 

where Ss,t *s the semi-group for the generalized flow associated with the velocity 
defined in ^ in the limit as io ^ oo. Furthermore the map lo —> JiJ* satisfies 
the invariance property 

SrT^Jx,y,t)^Srjx,y,t + T). (55) 
An immediate consequence of this theorem is 

Corollary 8 In the limit as £q oo, in the strongly and the weakly compress- 
ible regimes, as well as in the intermediate regime if the flow is non-degenerate, 
there exists a unigue invariant measure on Lf^^{M.'^ x ft) for the dynamics defined 

by^. 

The proof of Theorem |^ proceeds similarly as the proof of Theorem ||. In the 
non-degenerate cases, one study the convergence of (compare ^5|) 

e( r [ St',t+M^,s) - biy,s))dsY ^ 0, (56) 
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as Ti, T2 —00 for fixed x and t. The average in ^ can be expressed in terms 
of P, and it can be shown that the expression in ( ^6| ) converges as Ti, 
T2 — > —00 in the non-degenerate cases. In the strongly compressible regimes, 
because of the existence of a flow of maps, ^ is replaced by 

E{£\b{vU^),s)-b{vlM,s))dsy. (57) 

This average can again be expressed in terms of P, and it can be shown that the 
convergence of the time integral in |5^ depends on the rate at P looses mass at 
r = 0+ (i.e. the rate at which particles coalesce). The analysis of the equation 
in ^ shows that the process is fast enough in order that the integral over s in 
tends to zero as Ti, T2 — > —00 in the strongly compressible regime. In contrast. 



the equivalent of 57 in the intermediate regime in the e-limit can be shown to 
diverge as Ti, T2 —00. 

It can be shown that the invariant measure has finite correlation functions 
of all order, even though these results do not by themselves imply uniqueness 
of stationary solutions to the n-point Fokker-Planck equation. The task of 
studying the passive scalar is now changed to the study of the short distance 
behavior of these correlation functions. 



Incomplete self-similarity 

We finally turn to question Q4 and consider the scaling of the structure functions 
based on the invariant measure S-d* defined in Denote 

Sr.{\x-y\)=E\Sdt{x,y,t)r. (58) 

The dimensional parameters are Bq = B{0) ([temperature]^ [time] ^^), 
D ([length]2-«[time]-i), L ([length]). It follows that 

..„.(*P)"'";,.(I), 

where the /„ 's are dimensionless functions which cannot be obtained by dimen- 
sional arguments. For instance, the scalings in ^ correspond to different 
/„. It is however obvious from the equation ^ that, provided the limit exists 
and is non-zero 

lim 5„(r) = C„ (^^] = 0(r"(2-«)/2). (60) 

L— ►oo \ U J 

where Cn = liuir^oa fn{r / L) are numerical constants. The scaling in ^ is 
usually referred to as the normal scaling since, consistent with Kolmogorov's 
picture, it is independent of the forcing or the dissipation scales. In contrast, 
anomalous scaling is a statement that the structure functions diverge in the 
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limit of infinite forcing scale, L oo. In the spirit of Barenblatt-Chorin |Q, ||, 
we may say that normal scaling holds in case of complete self-similarity, whereas 
anomalous scaling is equivalent to incomplete self-similarity. 

It is interesting to discuss the existence or non-existence of the limit in ^ 
for both the coalescence and the non-degenerate cases. When the flow coalesces, 
because of the existence of a flow of maps and the absence of dissipative anomaly, 
the S2nS of even order 2n > 2 can be computed exactly It gives S2n{i') = co 
for n > C/(2 — ^), whereas 

S2n{r)=0{r<) for " < (61) 

where C, is given in Thus, for n < C/(2 — ^), 

/2„(r) = O ((r/L)«-"(2-«^ (g2) 

It follows that f2n and, hence, S2n tend to zero as i — > oo for 2 < n < C/(2 — ^), 
whereas they are infinite for all L for n > (2 — ^). In fact, in the coalescence 
case, it can be shown [|l8| that on scales much larger than the forcing scale L, 
the structure functions of order ji < C/(2 — ^) behave as 

52n(r) ~ C2„r"('-«) as r/L ^oo. (63) 

Thus in the coalescence case, it is more natural to consider the limit as i — s- 
of the structure functions, for which the expression in ^ shows the absence of 
intermittency corrections. 

In the non-degenerate case, one has 

52(r) =0(r2-«), (64) 

while perturbation analysis gives for the higher order structure functions ^, 

II 

52„(r) =0(r^="), (65) 
with C2n < n{2 - ^) for 2n > 2. It follows that /2(r) = 0(1), while 

/2„(r) = O ((r/L)^"-"(2-c)^ ^ 2n > 2. (66) 

In other words, as L ^ oo, 52 has a limit which exhibits normal scaling, whereas 
the iS2n's, 2n > 2, diverge. This may be closely related to the argument in |0, |[ 
that, in appropriate limits, intermittency corrections may disappear and higher 
than fourth order structure functions may not exist. We note, however, that 
Barenblatt and Chorin were discussing the case of infinite Reynolds number 
(here infinite Peclet number, k ^ 0) at finite L, whereas we require L oo. 
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